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Abstract 



We argue that the transverse mass spectra of identified hadrons as measured 
in gold-gold collisions at BNL's Relativistic Heavy-Ion Collider (RHIC) follows 
a generalized scaling law. Such a scaling behavior is motivated by the idea of 
a Color Glass Condensate, or more generally saturation of the gluon density. 
In particular, we describe the shapes of transverse mass spectra as a function 
of centrality. This scaling of the transverse mass spectrum is shown to be 
consistent with previously observed scaling of multiplicity with centrality. 
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I. INTRODUCTION 



In a previous work, two of us argued that the mean transverse momenta measured 
in relativistic heavy-ion collisions can be described by the intrinsic transverse momentum 
broadening seen at the Tevatron for pp collisions [I]. The mean transverse momenta of 
charged and identified hadrons increases universally with the square root of the multiplicity 
per unit transverse area for both pp and AA collisions. In addition, it has been shown that 
the total hadron multiplicity follows a scaling behavior motivated by the gluon saturation 
§• 

In this paper, we try to combine these two ideas in a consistent picture of transverse 
mass distributions of identified particles as a function of centrality. To do this, we apply 
theoretical ideas emerging in the Color Glass Condensate description of the high gluon 
density phase of QCD. Specifically, we describe the transverse momentum spectra using 
the recently measured data from BNL's Relativistic Heavy-Ion Collider (RHIC) for AuAu 
collisions at ^Js = 130 AGeV. 

The theoretical motivation for such scaling relations is as follows. At very high energies, 
the number of partons (primarily gluons) in a nucleus grows very rapidly. When the occu- 
pation number of these partons is large, they saturate [@-[5|] and form a novel state of matter 
called a Color Glass Condensate (CGC) P,|6HTT|j . The CGC has a bulk scale Q s (3> ^qcd) 
which is the typical intrinsic transverse momentum of the saturated gluons in the nucleus. 
The CGC can be probed in deeply inelastic scattering [[j3>[yl > * n photoproduction in periph- 
eral heavy- ion collisions [[14]] , in pA collisions | 15| , [16 l and in heavy-ion collisions P JT0|JTT[| . 

We remark, that scaling analysis of the transverse mass spectra at RHIC has been also 
pursued in fDJ although within a different physical picture. 



II. GENERALIZED M T -SCALING 
A. Scaling Relations 

Collisions of heavy-ions at high energies can be imagined as the collision of two sheets of 
colored glass, and a large multiplicity of gluons is produced. Since the occupation number of 
the gluons in each of the colliding nuclei is large, the problem can be treated classically fllP 



The only dimensionful quantities available to describe the collision are a) the transverse 
area of the two colliding nuclei a and b) the saturation scale Q 2 S we mentioned previously, 
which is determined by the density of partons in the transverse plane. This saturation scale 
is a function of both energy and centrality. All dimensionful quantities can therefore be 
expressed in terms of powers of one or the other scale times non-perturbative functions of 
the dimensionless product of the two scales. The initial momentum distribution of produced 



gluons in high-energy collisions is then found to be given by the relation [|TT 



dN g 



a drjd?pt 




(1) 



Here, f g is a universal, dimensionless function for the produced gluons which depends only 
on the ratio of the transverse momentum and the saturation momentum. Different energies 
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or scales are described by the same function with a correspondingly changed saturation 
momentum Q s . 

Motivated by the above relation, we want to test in this paper how well these scaling 
relations describe the actual data, i.e. the momentum distribution of produced hadrons. Ab 
initio, it is not obvious that there is such a relation as hadrons are produced nonperturba- 
tively at the deconfinement phase transition. There is also surely some transverse flow of 
the matter produced in these collisions and it is not clear to what degree this flow might 
distort the distributions produced early in the CGC. Also, hard processes are expected to 
modify the scaling relations at sufficiently high transverse momenta above the saturation 
scale Q s . 

We find that the transverse momentum distributions of identified particles are well de- 
scribed as a function only of the transverse mass m t = \Jrn\ + pf . This appears to be a good 
approximation so long as we are not too close to the mass of the particle in question. (At 
small values of p t we expect flow effects to be important which will invalidate this simple 
m t scaling form of the distributions.) We therefore parameterize the transverse momentum 
spectra as 

The universal function / depends on m t and the corresponding saturation scale for hadrons 
p s and incorporates m^-scaling. The constant Kh reflects the difference in abundances of 
various species of particles, and is of order one for the particles which we consider. We take 
it to be independent of p s . The hadron saturation momentum p s as well as a are parameters 
which are determined by the energy and centrality of the collision. We assume p s to have the 
same energy and centrality dependence as the saturation momentum Q s . Also, the function 
/ has to be extracted from the data. But once the universal function / is fixed for one 
data set, it should describe the other data by rescaling the parameters p s and a. Hence, the 
scaling relation implies that there is one and only one function which describes high-energy 
collisions at various centralities, energy and system size. This picture can certainly be valid 
only in some range of these parameters. It is the purpose of this investigation to determine 
how well scaling works for data from relativistic heavy-ion collisions. 



B. Comparison with RHIC Data 

We turn now to a more detailed discussion of the transverse momentum distribution 
focusing on identified hadron spectra as measured recently at RHIC in AuAu collisions with 
y/s = 130 AGeV |T8|-p3|. Note that in relation (0) all hadrons are described by only one 



function of the transverse mass m t = yrriQ+pf. The transverse mass is just the total 
transverse kinetic energy carried by that hadron. Plotting the transverse mass spectra for 
identified hadrons as a function of m t (not as a function of m t — mo) should then result in one 
single curve. This idea is not particularly new. It is similar to the m^-scaling behavior put 
forward by Hagedorn for pp collisions in his statistical model pl |. The universal function in 



the Hagedorn model is then a Bose-Einstein or Fermi-Dirac distribution and depends on one 
parameter, the slope parameter T s i ope . Contrary to this picture, we approach the problem 
in reverse and try to determine the universal function from the experimental data. 
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The transverse momentum spectra of pions, kaons, protons and antiprotons have been 



reported by the PHENIX collaboration at RHIC to quite large transverse momenta p t |19 
For the first time in hadronic collisions, one sees that the protons and antiprotons seem to 
be more abundant at higher transverse momentum p t > 2 GeV than pions. The slope of the 
transverse momentum spectra for protons and antiprotons is apparently much larger than 
that of the pions, so that the nucleon spectra overshoots the pion spectra at some pt > 2 
GeV. 

The scaling relation (0) can now be tested by plotting the transverse mass spectra. 
Figure [I] shows the preliminary minimum bias data of the PHENIX collaboration for charged 



pions, kaons, protons and antiprotons [I9[] and for neutral pions [5U| as a function of m t (not 
Pt). As one can see from Fig. [I], the data points for all hadrons seem to follow one curve over 
several orders of magnitude for a broad range of transverse mass. The neutral pion data is 
on top of the antiproton data within the error bars. But one notices also, that the curves for 
kaons and protons are shifted downwards and upwards, respectively. These are effects due 
to quantum numbers (strangeness and baryon number) which are not taken into account in 
the scaling relation. The production of strange particles is suppressed due to the strange 
quark mass, while the number of protons is enhanced due to the initial baryon number excess 
coming from the two colliding nuclei. Interestingly the proton-antiproton pair production 
seems to be not suppressed relative to pion production but kaon production is. We note that 
these changes are moderate, they shift the curves of the kaons down by about a factor two 
while shifting the proton curve up by about a factor two from the universal curve. These 
shifts correspond to the constant Kh in eq. |2|. Shifting the kaon and the proton spectra by 
these factors results in Fig. 0. Now all the data points are lying on one curve. It tells us that 
the shape of the spectra for the same value of m t are equal for pions, kaons, and nucleons. 
(Often times, a slope parameter is extracted to characterize spectra. This is extracted from 
the spectra close to the mass threshold, i.e. at m t = rrih- The larger the mass of the hadron 
nih the larger will be its slope parameter T s i ope , as the slope of the transverse mass spectra 
is increasing with mt- We can read from Fig. 0, that this effect will also continue for heavier 
particles but it is less pronounced as the slope of the curve levels off.) If scaling works for 
all particles, the transverse mass distribution of hyperons will then follow also the same 
universal m r curve of Fig. § modulo shape independent effects from baryon and strangeness 
number conservation which affect only the overall normalization of the distribution. 

Being more quantitative, the transverse mass spectra of Fig. can be fitted by a power 
law of the form (1 + m t /p s )~ n with the parameters n = 16.3 and p s = 2.71 GeV. We remark 
that the values for n and p s are strongly correlated, fits with different values for p s and 
properly adjusted n gives an equally good fit to the data points. The values for n and p s 
are constrained in such a way that they get about the same mean transverse mass of 

« = ^ (3) 

for each fit. Let us define now the local slope as 

1 d 



T s \ ope dm t 



\n(f(m t /p s )) . (4) 



Then we find that the local slope parameter for a power-law distribution in m t is given by 
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Aslope = — + -m t . (5) 
n n 

i.e. a constant term plus a term linear in the transverse mass. The constant is closely related 
to the mean transverse mass. The correction term proportional to the transverse mass 
originates from the non-exponential behavior of the power-law distribution which enhances 
the high-m t part of the distribution. With the fitted parameters given above, we find the 
following slope parameters for various hadrons at m t = m^. 175 MeV (tt), 196 MeV (K), 224 
MeV (p), 235 MeV (A), 247 MeV (H). The actual measured values for tt, K and nucleons are 
higher |TJ| as the experimental fits are taken around 0.5 GeV above the threshold. Hence, 



the apparently different slopes for hadrons can be also explained by a generalized m t -scaling 
of the transverse mass spectra. Note, that the relation (||) is similar to the one derived for 
radial flow of non-relativistic particles T slope = T Q + 0.5(3 2 -nth [f^J. Inpp collisions, the value 
for p s is smaller than in AA collisions, so that the deviations from an exponential shape of 
the transverse mass distribution are less pronounced. Therefore, one recovers in pp collisions 
the traditional m t -scaling behavior of Hagedorn's statistical model |24| with about similar 
slope parameters for all hadrons. We will discuss the centrality dependence below in more 
detail. 

We remark that the universal m t -scaling seen in Fig. Q does not rule out radial flow 
because at large momenta the hydro picture also predicts a universal curve for all hadrons. 
A clear indication of radial flow would be a deviation from the mt-scaling curve. These 
deviations should be especially pronounced for heavier hadrons. STAR measures at even 
lower p t than PHENIX and reports a larger slope parameter for antiprotons than PHENIX 
Z2J251]- This experimental finding cannot be explained by generalized m t -scaling. The 



flatter distribution at small m t close to compared to larger m t can be well described 
in hydrodynamical models with collective transverse flow [P5]-|2"£|j . In general, radial flow 
implies a violation of m t scaling, as the transverse momentum distribution then depends 
then on both, p t and m t . 

We expect that there should be radial flow effects in heavy ion collisions. These should 
show up in a relatively unambiguous way as the deviations from m t scaling for m t near 
threshold. What we try to describe in this paper are the gross features of the data for mt far 
from threshold. The consistency of the picture we present is simple within the Color Glass 
Condensate framework. It may turn out for reasons which we do not understand that radial 
flow preserves these simple patterns. For example, this might be plausible if the flow sets in 
very early in the collision. It might also turn out that the patterns we observe have some as 
yet undiscovered interpretation. The scaling relations we present nevertheless would remain 
as a simple phenomenological description of the data. 

Finally, we note that the universal mj-scaling for hadrons in heavy-ion collisions has 
been also observed at lower bombarding energies at CERN's SPS |29|j30|| but it has not been 
discussed in detail. In particular, no interpretation in terms of scaling has been given. 



III. SCALING WITH CENTRALITY 

The scaling relation (^) predicts more than the universal curve for the transverse mass 
spectra. As it should be valid for all centralities, the transverse mass spectra should scale for 
each centrality class. Moreover, there is only one scaling function and the transverse mass 



5 



spectra for different centralities can then be rescaled into each other by properly choosing 
the transverse area a and the momentum p s for each centrality bin. 

The first statement is probed in Fig. R] which shows the preliminary m r spectra of 



PHENIX for negatively charged pions and antiprotons for different centralities ||19|| . In- 
deed, it seems that there is a universal function for each centrality bin, even for the most 
peripheral one, which describes the pion and antiproton spectra simultaneously. The form 
of the curves looks similar for the various centralities, too. 

To check that the m t distribution is universal, we rescale the data points in absolute 
normalization and in transverse mass: 

1 dN h 11 dN h m t 

so that the data points for all centralities are lying on top of each other. We choose the 
most central bin as a reference curve. The scaling parameters A and A' are extracted for each 
centrality bin. The scaling parameter A is just u/a s {p s ) for a given centrality relative to the 
value for the most central bin. The parameter A' is the momentum p s for a given centrality 
bin divided by the one for the most central bin. Technically, we use the p^-spectra of charged 
particles as measured by the PHENIX collaboration |L8] to fix the scaling parameters A 
and A' for each centrality bin as the data has much better statistics. Then we use those 
parameters to rescale the mt distributions. 

Fig. |] shows the rescaled version of the transverse mass spectra for all centralities. The 
data points for all centralities are now lying on top of each other, i.e. scaling works reasonably 
well also for the centrality dependence of the m r spectra. Accordingly, also the charged 
particle p t spectra then can be rescaled for each centrality cut. There seem to be deviations 
from scaling at large m t in Fig. ^ as that the more central bins are suppressed relative to the 
peripheral ones. But within the error bars it is difficult to make a more definite conclusion. 
We see a similar trend at high p t for the rescaled charged particle spectra, but again within 
the error bars it is hard to draw a conclusion. 

One expects that a/a(p s ) rescales as a product of the transverse area times a logarithmic 
function of the centrality. Similarly, one expects p\ to scale as the cube root of the centrality. 
In the following we will check, if this is the case. 



A. Scaling of the Transverse Area 

The fitted scaling factor a/a s (p s ) for the different centrality bins are plotted in Fig. || 
versus the number of participants A part . We take the number of participants as reported 
by the PHENIX collaboration in [[31] for their centrality cuts. The scaling factors are given 



relative to the one for the most central bin at A part = 347. 

One sees from the figure, that the fitted scaling factors seems to have a more linear depen- 
dence on Ap art than to follow the curve for a transverse area A ~ A pa r t . Nevertheless, there 
is a correction factor coming from a s (p s ). The quantity a s is taken as a^ 1 ~ In (Ps/Aq CD ) 
and the saturation momentum is dependent on the number of participants A part |J. The 
factor a s induces then a logarithmic correction to the scaling factor. The scaling factor 
should then follow a curve proportional to the transverse area A divided by a s which is of 
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the form N pait In [p 2 / Aqcd) ■ For consistency, the argument of the logarithm is taken from 
the dependence of p s on Apart as will be discussed in the next section. The dependence of 
p 2 s with iV part has been fitted to be p 2 Jp\ c = 0.61 + 0.39(A r part /347) 1/3 (see figure |) and p s , c 
is the momentum scale for the most central bin. The dependence of a s with p s is then of 
the form a; 1 ~ In ((0.61 + 0.39(A^ part /347) 1/3 ) //i 2 ) where /1 2 is the ratio of A 2 QCD and p 2 c 
and adjusted to be 0.6. The corresponding curve is shown in Fig. || by the solid line and 
now follows very closely the extracted values for the scaling factor. 



B. Scaling of the Transverse Momentum 

The second scaling factor, p s , as extracted from the data is shown in Fig. || Again, 
the factors are normalized to the one for the most central bin. The centrality dependence 
of p s is quite weak for moderate to large number of participants but falls off for peripheral 
collisions. The fitted values of p s are compared with the expected behavior for the saturation 
momentum Q 2 S ~ A^t in the figure. It is seen, that p s changes less rapidly with centrality 

1 /6 

than Ap art . The reason is that the gluon densities at hand are still too small to reach this 

/ 1/3 X 1 / 2 

behavior fully. Fits of the form (c + Ap ai . t J describe the centrality dependence of p s much 
better (see figure and the previous discussion in section [111 A| ). Hence, as p s increases it 
appears to be tending towards the scaling with N^ t . 

Let us compare p s to the mean transverse momentum defined by 

, > P 3 s ■ C h /Ps d 2 x^Jx 2 - {m h /p s ) 2 f{x) D fm h \ 

(Pt) = 5—755 , 2 ,/ x = Ps ■ P\ — , 7 

Pi ■ L h /ps d 2 xf{x) V Ps J 

where is the vacuum mass of the identified hadron and P is some function of the ratio of 
the hadron mass and the scaling momentum. The momentum p s is proportional to the mean 
p t . In the limit Ap art — > 2 one should recover the mean p t and its corresponding momentum 
p s of pp collisions. The constant seen in the fit of p s with A^ part then stands for that finite 
mean p t (or p s ) already present in pp collisions. Note, that p s is not directly identical with 
the mean p t but differs by the factor P. As this factor P depends on the vacuum mass of 
the hadron, the mean p t for each hadron follows a different behavior with centrality which 
depends on the ratio of its vacuum mass to the momentum p s in a nontrivial way. 

The weak centrality dependence of p s is compatible with the small change of the mean p t 
with Apart seen at RHIC [ |T9"| , |3"2"] . The mean p t for charged particles measured in pp collisions 



at ^fs = 200 GeV has been measured by the UAl collaboration to be (p t ) = 392±3 MeV |33 
The energy dependence of the mean p t for pp collisions is known to be quite weak. STAR 
reports preliminarily that (p t ) = 508 ± 12 MeV for central AuAu collisions at y'i = 130 
GeV [p2| ,j3~^|. So, the ratio of mean p t squared for Ap art — > 2 to the value for central AuAu 
collisions, i.e. the value for p 2 at Ap art = 2, is calculated to be 0.60 which is in agreement 
with the results of Fig. ||. 

The centrality dependence for p s can be compared to the measured transverse energy per 
charged particle by the PHENIX collaboration |34] which is approximately proportional to 
the mean transverse momentum. The data points are flat from the most central bin down 
to Ap ar t ~ 76. The extracted values for p s in Fig. || follow closely this constant dependence 
on centrality for these values of Apart- 
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C. Multiplicity 



We can cross-check our scaling relations by computing the centrality dependence of the 
charged multiplicity. Integrating eq. (fj) over the transverse mass gives the multiplicity 

f = ^. Kh r d * xKx) = °A. KhF (^) , (8) 

dy a s Jm h / Ps a s \ p s J 

The dependence of the multiplicity on the scaling factors is analogous to the multiplicity 
dependence of initially produced gluons in the CGC (see e.g. and references therein): 

— (9) 
dy a s 

except for the last factor which takes into account normalization effects associated with 
different species of particles. These are small corrections to these formula, as e.g. the kaon 
spectra is lying about a factor two below the universal curve in m t (see our discussion of 
Fig. |j). Neglecting the dependence of the shift from rapidity y to pseudo-rapidity 77 on 
centrality, the charged multiplicity density reads 

1 dN ap 2 s 



iVpart dr] a s N pSLTt [ \ p s J \ p s J \ p s ) 

If we take the universal function / to be a power law in the transverse mass 



(10) 



then it can be analytically integrated to get 

f(^U ? I A Ps V~7l + (n-l)^ . (i 2 ) 

\p s ) (n-2)(n-l)\m h + p s J \ 1 'p.) 1 ' 

Indeed, one can perform a power law fit to the transverse mass spectra of the preliminary 
PHENIX data of the most central bin of Fig. |3] with n = 11.8 and p s = 1.65 GeV. We 
can now compute the charged multiplicity from the two scaling factors and the power law 
function for the most central bin. 

Quantitatively, the factor a/a s grows approximately like iV^t ln(c + iV,^) (see Fig. [|) 
and p 2 s grows like c + (see Fig. ||), so that finally we recover 

Npart dr] v " ) 

This form is the one expected from the high density QCD approach with a nearly constant 
contribution from soft scatterings and a term growing with In iVp art 0. It can be motivated 
by noting that at large iVp art it has a form consistent with our expectations concerning 
saturation. At low A^ part , the constant gets important reflecting non-perturbative physics, 
perhaps associated with a vacuum gluon condensate fl35| , |36]j . That the dependence of the 
charged multiplicity on iVp art comes out correctly is nontrivial. The multiplicity of particles 
is dominated by the contributions from small p t where it is possible to have deviations from 
the universal function due to resonance decays. As it turns out this does not seem to be the 
case qualitatively. 
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IV. TRANSVERSE MOMENTUM BROADENING REVISITED 



The eqs. (|7|) and (pi) give a relation between the mean transverse momentum and the 
charged multiplicity: 

(p t ) 2 ~ a s P 2 (m h /p s ) ■ ~ . (14) 

Except for the factor a s P 2 , this is the relation as studied in [[I]]. The mean momentum 
squared increases as the charged multiplicity per transverse area. This has been found to 
be in good agreement with the pp data from the Tevatron and the heavy-ion data from 
the SPS [|TJ. We expect a similar behavior at RHIC from scaling arguments. Compared to 
the previous work |J , we are now in the position to explain in more detail the effects of p t 
broadening, i.e. that the mean p t rises more rapidly for more massive particles compared to 
the purely kinematic behavior as seen in pp collisions at the Tevatron [37- §^] . The mass 



dependent factor P(m,h/p s ) gives a different steepness in the increase of mean p t squared for 
different hadrons, so that pions, kaons and nucleons have a different slope. As P{my l /p s ) 
increases with hadron mass m^, the slope will be larger for more massive particles. If the 
function / is of exponential form, / ~ exp(— m t /T), then the mean p t is given by 

W = T hr) K m (,n h IT) ■ (15) 

For heavy particles, <C T, one gets (pt) 2 ~ K jlm^T which is besides a small correction 
factor just the purely kinematic factor 2rrihT for non-relativistic particles. 

We show in Fig. ^| as a function of centrality the PHENIX data on the transverse energy 
per charged particle E T /N C EM. The data has large systematic errors, and taken alone 
would provide no evidence of p t broadening. The errors on this data nevertheless allow for 
consistency with the data on identified particle spectra. The STAR data on mean p t for 
charged particles shows a centrality dependence which appears consistent with our analysis 
22J32L 
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FIG. 1. Transverse mass spectra of identified hadrons for minimum bias gold-gold collisions as 



measured for y/s = 130 AGeV at RHIC (preliminary data taken from [IS, 2 
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FIG. 2. Transverse mass spectra of figure [l] where the proton and kaon/antikaon data points 
are multiplied by 1/2 and 2, respectively. 
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FIG. 4. Rescaled transverse mass spectra of ir and p of the previous figure. 
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FIG. 5. Scaling factor of the absolute normalization of the transverse mass spectra. All curves 
are normalized to the most central bin at A^ part = 347. 
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FIG. 6. Momentum scaling factor p s as extracted from the rescaling of the mt-spectra. The 
data for the transverse energy per charged particle from PHENIX [34] is plotted for comparison 
and normalized to the most central bin. 
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